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Abstract. INTEGRAL is the first gamma-ray astronomy mission with a sufficient sensitivity and angular resolution combina- 
tion appropriate to the detection and identification of considerable numbers of gamma-ray emitting sources. The large field of 
view (~ 30° zero response FWHM) enables INTEGRAL to survey the galactic plane on a regular (~weekly) basis as part of the 
core progr amme. The first s ource catalogue, based on the 1st year of core programme data (~5 Msec) has been completed and 
published ("Bird et al., 2004). It contained 123 y-ray sources (24 HMXB, 54 LMXB, 28 "unknown", plus 17 others) - sufficient 
numbers for a reasonable statistical analysis of their global properties. These were located to a positional accuracy of typically 
0.72 arc minutes. The detection of previously unknown y-ray emitting sources generally exhibiting high intrinsic absorption, 
which do not have readily identifiable counterparts at other wavelengths, is intriguing. The substantial fraction (roughly 20% 
of the total number) of unclassified y-ray sources suggests they must constitute a significant family of objects. In this paper we 
review the global characteristics of the known galactic sources as well as the unclassified objects with the twin aims of investi- 
gating how the unclassified set may fit into stellar evolution and improving our understanding of known X-ray binary systems 
through the non-thermal y-ray channel. In the context of the known systems we are very conscious that they constitute a y-ray 
selected set, and may exhibit subtle generic differences to the rest of the class. We present Log(N)-Log(S) distributions, angular 
distributions, and for systems with reliable distance estimates the spatial distributions within the Galaxy and luminosity func- 
tions. For the unknown sources, this statistical analysis has shown that they are most likely to be HMXBs containing a highly 
magnetised neutron star The lack of X-ray counterparts for these sources indicates a high degree of intrinsic obscuration. 

Key words. Gamma-rays: observations - X-rays: binaries ~ Galaxy: general - Galaxy: structure - Galaxy: stellar content 



1. Introduction and the First IBIS/ISGRI Catalogue 

INTEGRAL is the first gamma-ray astronomy mission with a 
sufficient sensitivity and angular resolution combination ap- 
propriate to the detection and identific ation of considerabl e 
numbers of gamma-ray emitting sources ( Wink ler et allEoO Si). 
With an observation time of roughly 1000 seconds required 
to detect a IMq neutron star emitting at the Eddington limit 
from a distance of lOkpc it is clear that a meaningful survey of 
discrete galactic gamma-ray emitting objects is possible in the 
lifetime of INTEGRAL. These will be located to a positional 
accuracy of typically 0.72 arc minutes. 

A significant fraction of the INTEGRAL core programme 
is devoted to regular scans of the galactic plane and a deep ex- 
posure of the galactic centre. This is facilitated by the large 
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field of view of the on-board telescop es (9° x 9 ° fully coded 
and 29° x 29° zero response for IBIS. IUbertini et al. (20^), 
which permits ^lO** second exposures of the galactic plane to 
be made on a regular basis (~every 12 days). The sample of 
sources used for this study is derived directly from the first 
IBIS/ISGRI catalogue llBird et alll2004l) . which represents the 
first year of the INTEGRAL Galactic Plane Survey (GPS) and 
Galactic Centre Deep Exposure (GCDE) observations. The 1" 
IBIS/ISGRI survey catalogue was compiled based on all Core 
Programme observations between revolutions 46 and 120 in- 
clusive. INTEGRAL/IBIS data is organised into short point- 
ings (science windows) of ~2000 seconds. OSA software ver- 
sion 3 was used to create images at science window level, in- 
cluding two iterations of processing to first produce a clean- 
ing catalogue, and a second processing to remove all ghosts 
from the images. The ~2500 individual science window images 
were then mosaiced using a custom tool to produce deep all-sky 
maps. An initial source list was obtained by searching all-sky 
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mosaics constructed in both 20-40 and 30-50 keV bands, us- 
ing the SExtractor tool. The 20-40 keV band was optimised for 
detection of soft sources, and a threshold of 7 sigma was ap- 
plied. The 30-50 keV band was judged to produce the cleanest 
images, and a threshold of 6 sigma was used. By inspection of 
the image statistics, it was determined that at maximum, one 
source in the catalogue would be a false detection. It has since 
been discovered that the indentification of IGR J 17460-3047 
was in fact a false detection. 

Positions for each source were provided using the centroid- 
ing function in SExtractor, which allows for a non-Gaussian 
PSF imposed by searching all-sky maps in the Aitoff pro- 
jection. The source positions were principally extracted from 
the 30-50 keV maps, which were least contaminated by back- 
ground structures. In the case of very soft sources, positions 
were taken from the 20-40 keV maps. Positions for previously 
known sources were cross-checked against archive positions, 
and found to match the expected point source location accu- 
racy welllGros et al., 2003). The typical point source location 
error range is 20"-3', with a 1 arc minute error circle for a 
source >10cr(90% confidence). The precision of the source lo- 
cations allow for the clear identification/association to sources 
seen by previous missions and hence source classification is 
possible. Additionally, those sources without previously seen 
high energy counterparts can be confidently identified as newly 
discovered INTEGRAL y-ray sources. Fluxes for each source 
were obtained by extraction of values from the all-sky mosaics 
at the positions of best fit provided by SExtractor. 

The 1st IBIS/ISGRI survey catalogue is constructed from 
approximately 5 Msec of observations spread throughout the 
period February 2003 to October 2003. The observations were 
clustered along the galactic plane; consequently most of the 
y-ray sources are stellar systems in one form or another. 
INTEGRAL observes the available galactic longitude in a sun- 
orthogonal manner dictating that most regions of the galactic 
plane have been observed intermittently during two typically 
three month periods within the annual cycle and spaced roughly 
six months apart. Since most of the galactic y-ray sources are 
highly variable the catalogue is a time averaged map and the 
flux value quoted depends on its total accumulated flux aver- 
aged over the period of the observation for that region of the 
sky, no other attempt has been made to correct the recorded 
flux to an "average" value. A number of transient objects are 
therefore missed if their cumulative emission over the survey 
period does not rise above the aggregate threshold. 

Of the 123 sources in the catalogue, 5 were unambiguously 
identified with AGN, 5 with white dwarfs, 4 with radio pul- 
sars (free neutron stars), 3 supernovae remnants and one clus- 
ter Of the remainder, 54 were identified with known low mass 
X-ray binary (LMXB) systems, 24 with high mass X-ray bi- 
narys (HMXB) and the remaining 28 sources have no firm 
classification and are here after described as "unclassified". In 
this paper we have only considered the LMXBs, the HMXBs 
and the unclassified objects, discounting IGR J17460-3047 as 
a flase detection; the small numbers of the remaining object 
types preclude any meaningful statistical analysis. Formally, 
the LMXBs w ere classified as having the primary star with a 
mass < I Mq jLiu et al.ll200ll) and HMXBs having a primary 



star with M> 1 M© (Liu et allboOll) . Since the distribution of 
the unclassified sources is strongly associated with the galactic 
structure we have assumed, for the purposes of deriving LogN- 
LogS distributions that the contamination of the sample by ex- 
tragalactic sources is extremely small, or even non existent. 

The number of galactic gamma-ray sources in the first 
IBIS-ISGRI catalogue are compa rable to the early X-ray cat- 
alogues e.g. IWarwick etai1(ll98ll) and are sufficient for a pre- 
liminary investigation of their global characteristics on a statis- 
tical basis. It should be noted that the sources presented in the 
IBIS-ISGRI catalogue are gamma-ray selected and, although 
many of them are well known X-ray objects, it is not surpris- 
ing that there are considerable differences to catalogues of X- 
ray sources of corresponding depth. This is due to the fact that 
IBIS-ISGRI detects photons at energies largely unaffected by 
photoelectric absorption in universal abundance matter Hence 
objects highly photo absorbed in the classical X-ray band are 
clearly visible at gamma-ray energies, a factor confirmed by the 
significant fraction of sources previously unclassified. There is 
clearly a thrust to find out what these unclassified objects actu- 
ally are and how they fit into the overall picture of stellar evolu- 
tion. There are two main methods to achieve this end. One is to 
perform follow-up observations on the individual gamma-ray 
sources at other wavelengths in order to identify a counterpart 
and thus permit more detailed studies, where possible, over a 
wider spectral range. The other approach is to study the sample 
on a statistical basis and compare their overall characteristics 
with other known generic sets of astronomical objects. There 
is a strong degree of complementarity in the two methods; here 
we investigate the global characteristics of the y-ray sample. 

2. The Log(N)-Log(S) Distributions 

In this section we construct the number flux distribution for 
the various source types. We have assumed the usual power- 
law form of N(>S) =KS " for the relationship. Figure^shows 
the raw form of these relationships for all galactic sources, the 
HMXB, the LMXB and the unclassified sources separately. 
The source intensities have been converted into 20-40 keV 
Crab equivalent units. The decrease in source density below 
a few milliCrab is due to the flattening of the absolute source 
distribution as the sample is flux limited around that value. This 
corresponds to an energy flux of ~3 x 10"" erg s"' cm"^ for a 
Crab-like spectrum. 

As the sky coverage of the first survey is not complete, 
these curves are only indicative of the true galactic distributions 
and must be corrected for both the area of sky covered and the 
depth of exposure (or the minimum detectable flux (MDF)) at 
each point. This is not straightforward due to the residual sys- 
tematic structures in the sky maps, however it is possible by in- 
spection of the error and exposure distributions to parameterise 
the variation in MDF as a function of exposure, thereby taking 
into account the effects of the residual systematic variations. 

From this relationship and the general sky exposure map 
we can then construct Figure |2l which shows the area of sky 
(in square degrees) observed as a function of MDF (6 sigma in 
Crab units). 
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Table 1. Slope of Number/Flux relation 



All galftctic sources 



LMXBs 




0,001 



0,010 0,100 
Flux (Crab) 



1,000 



Fig. 1. The raw Number-Flux relationship for the sources in 
the 1st INTEGRAL/IBIS survey. The curves shown are for. 
All Galactic Sources; LMXB's, HMXB's and the unclassified 
sources. 

We can then use this curve to correct the observed logN- 
logS resulting in the relationships shown in Figure|3] 

The limiting detection threshold of 6 sigma is suffi- 
ciently high so that the ma ximum likelihood (ML) method 
llMurdoch & Crawford! 197? ) can be used to calculate the best- 
fit values of the slope of the number-flux relationship without 
the uncertainty in the correction factor for weaker sources dom- 
inating the correction itself. This has been performed for both 
the uncorrected and corrected distributions and the value of the 
slope found in each case is shown in Table 1 , 

It is of interest to consider and compare the Log(N)-Log(S) 
distributions of the various sub-groupings. The -0,91 slope of 
the power law for all the combined galactic sources is close to 
the -1 value expected for a uniform infinite plane distribution 
and lies between the value (-0,79 + 0,07) m easured for galac- 
tic sources by ASCA (Sugizaki et allEoOll) in the classic 2-10 



Source Class 


Uncorrected Slope 


Corrected Slope 


All Galactic Sources 


-0,81±0,10 


-0.91+0.09 


LMXB 


-0,81±0.13 


-0.95±0.13 


HMXB 


-0,65±0,15 


-0.81+0.15 


Unclassified Sources 


-1.79±0.37 


-2.11+0.46 



keV X-ray band and th e (-1.1) slope der ived from the Einstein 
galactic plane survey ( iHertz & Grindlav„ .1984.1 , It is signifi- 
cantly steeper than the value ( — 0,5) derive d for the bright 
galactic sources from the Ariel 5 catalogue dWarwick et all 
.1981.) , This all-source value is clearly dominated by the numer- 
ically superior LMXBs, which when detached provide a slope 
of -0,95 ± 0, 13 as a sub-set. The value for LMXBs is thus con- 
sistent with - 1 , as qualitatively expected for a population with a 
larger scale height ( Sugizaki et al,, 200 ll) and less aligned with 
the spiral arm tubes. In contrast HMXBs are typically located 
within spiral arms and this would automatically lead to a value 
of a closer to 0,5, In this context, the -0.81 +0.13 slope mea- 
sured for the HMXBs appears slightly flatter than the LMXBs 
and may reflect the likely location of these objects along the 
spiral arms, however the precision of the LMXB and HMXB 
slopes is insufficient to distinguish between the two groups. 
The steep slope (-2.11) estimated for the unclassified sources 
is intriguing and could imply a dilferent source population at 
low fluxes, this is far steeper than we may expect from even an 
extragalactic contamination. However, it must be remembered 
that there are relatively very few of these objects, and they are 
also amongst the weaker sources in the survey therefore care 
must be taken in any interpretation of this slope, as indicated by 
the large uncertainty obtained. Given the angular distribution of 
the unclassified sample (see Figure |6j, which indicates a con- 
centration towards the central regions of the galaxy it could be 
argued that this steeper distribution may be due to an increase 
in their space density towards the galactic centre. 
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Fig. 2. The area of sky observed in the survey as a function of 
file minimum detectable flux (6sigma), The flux is in equivalent Fig- 3. The Number-Flux relationships as shown in Figure 1 but 



Crab units. 



corrected for exposure and sky area observed. 
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Fig. 5. The galactic distribution of sources in the 1" IBIS/ISGRI catalogue. Filled circles represent HMXBs; open circles rep- 
resent LMXBs; triangles represent unclassified INTEGRAL sources. Overlaid are contours of exposure time, 50, 200 and 700 
ksec. 



If we follow this line of reasoning, and if we assume that 
the newly discovered INTEGRAL unclassified sources are all 
located in the galactic bulge region we can investigate their spa- 
tial distribution by assuming the distribution is the same in the 
line of sight towards the galactic centre as it is orthogonal to our 
line of sight. We set the zero point to lie at the peak of the distri- 
bution and take absolute values of the distance of sources from 
this point because of the low number of statistics. The resulting 
distribution is shown in Figure 0] The decay in the number of 
sources away from the Galactic Centre can be fitted by a power 
law with an index of ~L3 + 0.8. This distribution, coupled with 
the fact that the sensitivity of the first catalogue does not extend 
appreciably beyond the galactic centre region, would imply that 
for a disk population the sources would obey a LogN-LogS dis- 
tribution with a slope of 1.7 + 0.4, only about 1 standard de- 
viation away from the observed value quoted above. In general 
due to the limited number of sources in the first catalogue, and 
also the remaining systematics affecting the imaging process 
these results are necessarily subject to several sources of error 
which will be overcome with the release of the second cata- 
logue. Nevertheless they provide qualitative information on the 
behavious of galactic source populations above 20 keV. 



3. The Angular Distributions 

The distribution of the y-ray emitting sources on the sky is 
shown in Figure 13 The open circles correspond to the posi- 
tions of LMXBs, the filled circles to HMXBs, the unclassified 
sources are represented by the triangles. Superimposed on this 
sky map are contours of the sky exposure. It is readily appar- 
ent that the first year of INTEGRAL observations have been 
strongly biased towards the study of the Galactic Plane and 
in particular the Galactic Centre, hence the first IBIS/ISGRI 
catalogue is essentially a register of galactic sources. Figure|5] 
also documents the lack of uniformity of the exposure along 
the Galactic Plane. Despite this non-uniformity, the concentra- 
tion of the y-ray selected HMXBs along the Galactic Plane is 
clearly apparent, and likewise the grouping of LMXBs in the 
galactic Bulge region. 

Figure|6lreveals the considerable differences in the longitu- 
dinal distribution of the HMXB and LMXB sources. Whereas 
the LMXB are concentrated in the Galactic Centre region, 
the HMXB are spread more extensively along the Galactic 
Plane, but in a non-uniform manner, with some evidence for 
bunching in the vicinity of the tangential points of the spiral 
arm struct ure. The locations of th e spiral arm tangents were 
taken from lEngknaier et alJ (Il999l) . Additionally, a number of 
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HMXBs are seen at l" ~ ± 90", which is indicative of sys- 
tems located within our particular spiral arm. This association 
had been previously noted through inspection of the Ginga 
data bylKo vama et al. ( 1990) and in the RXTE/ASM data by 
lGrimmetal] ( i2002l) . 

The association with the spiral arms is entirely to be ex- 
pected since the high mass binaries are young stellar systems 
and should be attached to regions where star formation has re- 
cently taken place, such as the spiral arms. We will come back 
to this point later, when we discuss the spatial distribution of 
the y-ray HMXBs. 

Upon first inspection the unclassified INTEGRAL sources 
appear to have a longitude distribution concentrated around the 
Galactic Centre, similar to that of the LMXBs. To some extent 
this is misleading as the Galactic Centre has had the most sky 
exposure and hence we are more likely to detect new systems 
in this region. However upon closer examination, whilst some 
emulate the distribution of LMXB, it can be seen that the un- 
classified sources do not have a symmetrical distribution about 
the Galactic Centre and show some tendency to cluster around 
the locations of spiral arm tangents, specifically the Scutum, 3- 
kpc and Norma arms, as shown in Figure|6l Specifically the un- 
classified sources appear to concentrate in the 3-kpc and Norma 
spiral arms. 

The angular distribution off the Galactic Plane of the 
sources is seen in Figure The HMXBs exhibit a much nar- 
rower range of angular separations as expected if they are con- 
fined to the Galactic Disc. The much wider spread seen in the 
LMXBs is indicative of a population derived from the Galactic 
Bulge. Predictably the higher spread in the latitude distribu- 
tion of the LMXB as opposed to the accumulation of HMXB 
closer to the Galactic Plane is to be expected from the relative 
ages of their progenitor stars. The unclassified sources appear 
to have a sharp latitude cut-off similar to that of the HMXBs, 
this may indicate that they too are a primarily a population from 
the Galactic Disc. The unclassified sources, however, do not 
precisely conform to either of the subgroups; this is not unex- 
pected as they could consist of a mixture of both, or be a sep- 




Number of sources 



Fig. 4. A histogram of the absolute distance of unclassified 
sources away from the galactic centre perpendicular to the line 
of sight (assuming they are ~8 kpc away). The dashed fine rep- 
resents a power-law fit to the histogram. 
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Fig. 6. The Galactic longitude distribution of HMXBs, LMXBs 
and unclassified sources. Labelled are the locations of the spiral 
arm tangents: (a) - Sagittarius; (b) - Scutum; (c) - 3-kpc; (d) - 
Norma; (e) - Centaurus. 

arated generic set. Although the longitude distribution is not 
entirely dissimilar to the LMXBs, this is probably distorted by 
an effect of exposure and the source concentration around spiral 
arm tangents suggests an association with high mass systems. 
This association is further supported by their latitude profile 
that more closely resembles that of the HMXB. However the 
unclassified sources are very much associated with the Galactic 
Centre region, this association may be generic in some way and 
possibly an artefact of the selection of the sources through the 
y-ray channel, which is not sensitive to the effects of higher 
photoabsorption generally suffered through the conventional 
X-ray observations. 



4. The Spatial Distributions 

The distances to a number of the HMXBs and LMXBs 
were found through an examination of the availabl e data 
as published in the litera ture ( White & van Paradiii, 1 19961; 
iMacomb & Gehrelslll99 9': Kuulk ers et all i2003). thus allow- 
ing the construction of a 3D model of the source distribution. 
An updated list of distances collected from available litera- 
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Fig. 8. The galactic distribution of the subset of HMXBs (filled circles) and LMXBs (open circles) for which distance estimates 
are available in the literature. The Sun is represented by the star symbol and is 8 kpc away from the galactic centre. Superimposed 
is the 4-arm sprial model of .Vallea(i2002.) . Source distance errors are on average of the order of a kpc. 
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Fig. 7. Angular distribution of sources away from the Galactic 
Plane. Note that the number of LMXBs is divided by 2. 

tur e for neutron stars and bla ck holes in LMXB as reported 
by Ijonker & NelemansI (l2004l) represented one key source of 



distance information. Different methods were used in deriving 
the distances of such objects. For neutron stars characterized 
by the presence of type I bursts, distance can be inferred from 
observations of Eddington limited bursts. Type I bursts in Low 
Mass X-ray Binaries are thermonuclear flashes on the neutron 
star surface jLewin et aUil99 3). During some type I bursts the 
energy release is high enough such that the luminosity at the 
surface of the neutron star reaches the Eddington limit (for a 
typical neutron star the Eddington luminosity value ranges be- 
tween 2 - 4 X lO-'^ ergs cm"^ s"' ). During these events the 
neutron star's atmosphere expands because of radiation pres- 
sure. In this expansion and subsequent contraction the lumi- 
nosity stay constant near the limit and assuming a pure black 
body radiation there is a drop in temperature during the ex- 
pansion of photosphere radius. The enormous number of bursts 
and new burst sources discovered with the WFC s/BeppoSAX 
gave a ver y powerfu l tool in t his respect (see e.g. 'Coc chi et alJ 
(l2001„) : .Kuulkers et al.l2003.) : lin't Zand et al..(.2004.) ). In other 
cases the spectral pr operties of the companion star have been 
used. Where possible Jonker & Nelemans ( 2004) compared the 
two diff'erent methods and the discrepancies seem to indicate 
that a larger distance is derived from the first method while us- 
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ing the second approach the distances are underestimated, pos- 
sibly indicating erroneous spectral classification of the com- 
panion star due to fast rotation. For the black hole candidates 
in our list three different methods have been used to derive dis- 
tances, namely by using the interstellar absorption properties, 
from the observed proper motions of receding and approaching 
lobes when symmetric jets are present, and from the compar- 
ison of the absolute magnitude of the donor star with the ap- 
parent magnitude using radius, spectral type and luminosity as 
determined directly from observation al data . 

A model of the galaxy taken from i Wui3 ( l2002h is shown in 
Figure|8]with the location of HMXB (filled circles) and LMXB 
(open circles). It is immediately obvious that the HMXBs and 
LMXBs do not share the same spatial distribution. The LMXBs 
can be seen to cluster at the Galactic Centre as implied by 
their angular distribution discussed in the previous section. 
The LMXBs are clearly clustered in the central ~3 kpc of 
the galaxy, while the HMXBs are found outside that radius, 
in the region associated with the spiral structure. However, 
for those systems >3 kpc from the Galactic Centre, it is im- 
possible to say whether LMXB or HMXB follow the spiral 
structure more strongly. Theoretically the appearance of the 
HMXB epoch should lag behind the leading edge of spiral 
arms. This is brought about by the density wave inducing star 
formation and the HMXB phase subsequently developing some 
time later. A broad time window is to be expected since the pe- 
riod of elapsed time from star formation to the appearance of a 
HMXB system is dependent upon the initial masses of its com- 
ponents, and would naturally re flect this spread. For example 
iTauris & van den H euveJ ('2 0041) calculate an age of 24.6 mil- 
lion years for a sys tem w hich start s with 14.4 Mq and 8 Mq 
stars, and Ivan den Heuvel & Heisd ( 1 19721) calculate an age of 
15 million years for a system which starts with a 16 Mq and 3 
M0 star 

An investigation of this systematic time lag is beyond the 
scope of this paper. An analysis of this would require a larger 
population of HMXBs with well defined distance measures and 
a model of the galaxy which incorporat es the qua si-stationary 
spiral structure hypothesis proposed bv iLin et alJ lfl969) . The 
key parameter in modelling the spiral structure is Q.p the an- 
gular velocity of the spiral pattern. Recent measu rements of 
Qp ba sed upon observations of open clusters by Dias & Lepind 
( I2OO5I) strongly indicate that the spiral pattern rotates as a rigid 
body with Qp ~25 km s ' kpc 'at the distance of the Sun from 
the galactic centre. 

The scale heights of the HMXB and LMXB systems are 
shown in Figure and Figure |5J). In each case a simple ex- 
ponential decrease in numbers as a function of the absolute 
distance h "above" the galactic plane, N = k e""'', was used 
to model the data sets. Making the simplistic assumption that 
all the unclassified sources reside close to the galactic bulge 
region, we constructed a pseudo scale height distribution for 
these unclassified objects as shown in Figure|9l;. Table|2]shows 
the associated values for the scale height distributions of the 
three catagories. 

It is anticipated that the scale heights of the HMXB 
and LMXB populations are different, and that the vertical 
distributions reflect the age disparity of their parent stars. 



a) HMXB 



z ri«PDi 



b) LMXB 
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Fig. 9. The distribution of the height off' the plane of: a) 
HMXBs; b) LMXBs; c) Unclassified sources. The dashed lines 
represent exponential fits. 
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Table 2. Summary of the chracteristics of the principal popula- 
tions of the IBIS-ISGRI sources in terms of their distributions 
above the Galactic Plane. 





k a Scale Height (pc) 


HMXBs 
LMXBs 

Unclassified sources 


9 ± 3 4 ± 1 240!^^ 
13 ± 3 2.9 ± 0.6 350^_ll 
12 + 3 4 ± 1 240!^° 



we adjusted r^, to match our measured scale heights. All dis- 
tances are in kiloparsecs. The standard galaxy model uses a 
mass ratio of 2:1 for disk:bulge. To correct the LMXB popu- 
lation both model components were used whereas the HMXB 
population uses only the disk component. 

To correct the observed luminosity function we calculate 
the fraction of the galactic mass which is observable at a given 
luminosity. The correction is taken from iGrimm et alJ toO^ 
and is shown below: 



Unfortunately, the precision of the observed scale heights 
shown in Table|2is such that the HMXB, LMXB and unclas- 
sified sources have scale heights which are within one sigma 
of each other and are hence indistinguishable. However it is 
interesting to note that the scale height of this y-ray selected 
HMXB set is consistent with the value of 150pc found for th e 
HMXB selected from RXTE/ASM data (' Grimm et all 12002). 
The lack of errors on the HMXB scale height as measured by 
iGrimm et al. ( 2002) make it difficult to accurately compare the 
values. Likewise the vertical distribution of the y-ray selected 
LMXBs is compatible with the 410+^™' pc scale height of the 
RXTE/ASM sample of the LMXB disk population however, 
our HMXB scale height lies only ~2cr from this measurement 
as well. If we take the unclassified source "scale height" at face 
value, then it is very interesting to note that, whilst they ap- 
pear to conglomerate in the galactic bulge region, their vertical 
distribution is the same as the HMXB systems, reflecting the 
similarities found in their longitude distributions. 



5. Luminosity Functions 

Many the sources in the first IBIS/ISGRI catalogue have good 
distance measurements and it is possible to construct their asso- 
ciated luminosity functions. However corrections are required 
to account for non-uniform exposure, the flux limited nature 
of the sample and any incompleteness in identifying object 
classes. This correction is performed by constructing a model 
of the stellar mass d istribution of the galaxy as performed by 
jGrimm e t al. ','20021). A galaxy model consisting of a disk and 
bulge component was constructed using the bulge and disk 
models of lDehnen et alJ lll998l) : 



PBiilge - PO,Bulge ' 



PDisk - PO,Disk ' exp 



r 

rd 



■ exp 



r2 + 



7^ \ 



(1) 



(2) 



Where po,Buige and po.oisk are the normalisations, r is the 
distance in the plane to the galactic centre, z is the distance out 
of the galactic plane and R is the distance from the galactic 
centre in spherical coordinates. The other parameters are: q is 
the oblateness of the bulge; ro is the scale length of the bulge; 
r, is the truncation radius of the bulge; r^ is the scale length of 
the disk; i, is the vertical scale of the disk; r,„ is the inner disk 
cut-off. These parameters were taken from fit of RX TE/ASM 
XRBs to the galactic model bv lGrimm etaP (l2002l) . However 



dN 
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dLj^^^^ M{< D{L)) 



(3) 



Where ^ is the true luminosity function, is the ob- 

served luminosity function, M,o, is the total mass of^the galaxy, 
M(<D) is the mass of the galaxy inside a distance D from the 
Sun. D(L) is defined as: 



D{L) — min 



4nF, 



-,D„ 



lim 



(4) 



Where F;„„ is the limiting flux of the sample and D,„a v is 
the maximum distance from the Sun of the sources used in 
constructing the luminosity function. We take F/,„, as ~5mCrab 
in the 20-100 keV energy band which corresponds to ~8.5 x 
10"" ergs cm"^ s"' . The maximum distance at which we detect 
sources approximates to the galactic centre distance. As a re- 
sult this means that it is only at relatively low luminosities that 
DflMj: does not correspond to D(L). As previously discussed the 
majority of our exposure is in the vicinity of the galactic cen- 
tre. To compensate for the uneven exposure we modified Fum 
for the four quadrants of the galactic plane. As the sensitivity 
is proportional to the square root of the exposure and that the 
maximum depth of our sample is ~8 kpc the change in the lu- 
minosity function correction within our range of exposures is 
<10% 

Accreting X-ray binaries are the brightest class of X-ray 
sources and their emission properties are determined by the 
nature of the compact object i.e. BHC or neutron star The 
strength and geometry of its magnetic field on the neutron stars 
are also important as is the geometry of the accretion flow from 
the companion. A neutron star with a strong magnetic field 
(~10^^ Gauss) will in fact disrupt the accretion flow at sev- 
eral hundred neutron star radii and material is channelled to 
the magnetic pole. In the case of low magnetic field (< 10'" 
Gauss) the disk may approach close to the neutron star sur- 
face. HMXBs typically have the accreting compact star as a 
high magnetic field NS or BH and LMXBs have typically a 
low magnetic field NS or BH. Strongly or weakly magnetized 
neutron stars are identified by the presence or absence of pulsa- 
tions in the light curve of the system. X-ray pulsars are seldom 
in LMXB, which generally show X-ray bursts that are sup- 
pressed when B is > 10" Gauss jLewin &.TossL 1 19831) 1 and 
this also explains why bursts are not detected in HMXB sys- 
tems. For a recent review on acc reting neutron star and black 
hole properties see |Psaltislil2004 . which also lists the energies 
of cyclotron lines in accretion powered pulsars used to deter- 
mine the values for the magnetic field. Magnetic field values 
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O Unclassificd- 




aj Black hole 



Fig. 10. The 20 - 100 keV luminosity functions of: a) LMXBs; 
b) HMXBs; c) Unclassified sources. The thin line represents 
the uncorrected luminosty function. The thick line represents 
the corrected luminosity function as described in the text. The 
unclassified sources are uncorrected as their source population 
is undetermined. 



for X-ray pulsars in the RXTE/BeppoSAX era are reported in 
ICoburnetar ('2002V 

Figure [Tol shows the cumulative y-ray (20 - 100 keV) lu- 
minosity functions for LMXB, HMXB corrected as outlined 
in the above text. Also included is a "luminosity function" for 
the unclassified sources, assuming that they are all at a com- 
mon distance of 8 kpc. It is noticeable that, whereas the low 
mass systems spread their 20-100 keV luminosities relatively 
consistently up to and beyond 10^^ ergs cm"^ s"', the major- 
ity of the high mass systems have a strong inclination to reside 
around lO^*" ergs cm"^ s"'. The distribution of the unclassified 
sources exhibits the low luminosity characteristics of the low 
mass systems, but also a tendency to favour luminosities close 



b) Neutron star 




Luminosity (ergs/s) 

Fig. 11. The 20 - 100 keV luminosity functions of the 
IBlS/ISGRl sources separated according to the nature of the 
compact object: a) black holes; b) neutron stars. The thin line 
represents the uncorrected luminosty function. The thick line 
represents the corrected luminosity function as described in the 
text. 



to 10^^ ergs cm"^ s"' . Whereas FigurefTOlpresents the y-ray lu- 
minosity functions of the known objects selected by the donor 
system. Figure ^2 shows the luminosity functions of known 
sources grouped by the nature of the compact object, black 
hole systems and neutron stars. There are marked differences 
in their individual distributions. Whereas most of the 20 - 100 
keV luminosities of black hole systems exist around lO-'^ ergs 
cm"^ s"', neutron stars do not emit y-ray fluxes much in excess 
of 10^^ ergs cm"^ s"'. Although subdividing the classifications 
has resulted in dwindling numbers of examples and hence sta- 
tistical significance, it is interesting to note from Fig.llOt and 
Fig. that the unclassified sources are most likely neutron 
star systems with a significant fraction, but probably not afl, 
cohabiting with high mass companions. The use of a ~10^^ 
ergs cm"^ s"' soft y-ray boundary to separate neutron st ar and 
black hole systems is not a new concept. [Barret et al] (Il996l) 
found for high luminosity-high energy emitters, the ~100 keV 
luminosity of black hole systems is commonly in the range 1 - 
6 X lO-'' ergs cm"^ s"', whereas for neutron star systems it is 
<10^^ ergs cm"^ s as confirmed by INTEGRAL 



6. Summary and Discussion 

The INTEGRAL galactic plane survey and ensuing first cat- 
alogue has offered for the first time the possibility to inves- 
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tigate the global y-ray characteristics of galactic sources on 
a reasonable statistical basis. In this paper we have studied 
the generalised observational features of y-ray selected previ- 
ously unclassified objects, and compared them with the cor- 
responding parameters of known HMXB and LMXB systems. 
Although the sample of sources has been y-ray selected the 
Log(N)-Log(S), angular and spatial distributions of the known 
classes, as expected , emulate those of X-ray selected sets e.g. 
iGrimm et al.l ( l2002l) . there appear to be subtle differences. Of 
great interest however are the 27 "unidentified" objects, which 
naturally have overhanging questions such as what are they and 
how do they fit into mainstream stellar evolution? 

Clearly a vigorous programme of follow-up studies on 
individual objects, using such missions as XMM-Newton, 
Chandra, and RXTE for timing studies, is one means to 
pursue the answers t o these ques tions, a nd such activitie s 
are already underway . [ Walter et alJ Eo03 ) ; Hi ll et alJ ilOO^: 
[Steohen et al.' (EoO?); 'Revnivtsev' (^2003); M att & Guai nazzi* 
(2003); Rodriguez et al. (2003), and a number of identifications 
with objects that emit at other wavelengths have been made. 
The complementary statistical approach employed here essen- 
tially comes to the same conclusion i.e. that a large fraction 
of the unclassified objects are obscured high mass X-ray bi- 
nary systems. Detailed studies of the main distributional char- 
acteristics of the unclassified sources when compared to the 
equivalent distributions of known X-ray binaries all point in 
this direction. Their Log(N)-Log(S) and broad angular distri- 
butions are consistent with them inhabiting the central 3-4 kpc 
region of the galaxy, similar to many of the LMXBs. However 
a closer look at the angular profile shows that it is not the cen- 
tral bulge that determines their finer scale longitude distribution 
but instead the structure of the inner spiral arms, indicating 
a considerable number of the unclassified objects are associ- 
ated with younger stellar systems, and hence probably HMXB. 
Likewise the galactic latitude distribution of the unclassified 
sources strongly reflects that of HMXB. More tenuously, the 
luminosity function of the unclassified sources, assuming them 
to be close to the galactic centre, exhibits features that are sim- 
ilar to HMXB systems and possibly those incorporating highly 
magnetised pulsars rather than black holes as the compact com- 
panion. 

It is hardly surprising that INTEGRAL should discover 
a population of previously unnoticed sources. The catalogue 
is y-ray selected, and INTEGRAL operates above the energy 
threshold for which significant photoabsorption takes place 
in universal abundance matter, so that a strong emitter above 
~30keV can be rendered insigni ficantly weak in the c lassi- 
cal X-ray band. As discussed by iLutovinov et alJ (12005:) . the 
HMXB systems constitute the most likely candidates. The stel- 
lar wind accretion mechanism that dominates in HMXB sys- 
tems, as opposed to the Roche lobe overflow associated with 
LMXBs, naturally provides a suitable dense and strongly ab- 
sorbing circumstellar wind to veil the X-ray emission in a sim- 



ilar manner to the ( 


:ase of the Seyfert 2 AGN configuration 


dTumer et al.l 


■boontl 


lMali7iaetal.[l200nh. 



Exactly how the y-ray selected highly absorbed systems fit 
into the overall picture of binary star evolution is currently un- 
clear. This is principally due to the requirement of a suitable 



companion able to provide a sufficiently strong stellar wind 
mass loss rate (M > 10""^ M q yr"') that is capable of gen- 
erating a suitably dense gas surrounding the compact object to 
stifle the X-ray emission. We need to understand why they are 
diff'erent to " normal" HMXB. Do these highly absorbed sys- 
tems relate to the mass/giant nature of the primary star, or to the 
orbit configuration, or are they experiencing a phase the binary 
systems routinely pass through; but have not been exposed at 
other wavelengths? Do they need to be HMXB? Clearly a series 
of dedicated observations are required to solve this problem. If 
they are not all HMXB systems, then could some of them be 
intermediate mass X-ray binaries (IMXBs) dvan den Heuvell 
J975). Some of the sources could be systems t hat are passing 
through, or close to the co r nmon envelope phase Jlben & Liviot 
ll993UTaam& Sandauisl l2000l) . However statistical grounds 
this seems unlikely due to the requirement that their relative 
numbers need to compatible with the relatively short timescale 
(< 10^ yr) of this phase, which is dictated by the drag force 
created as the compact star moves through the envelope of the 
extended companion. 

Clarification of the above general issues raised by the 
INTEGRAL source discoveries will also facilitate a better un- 
derstanding of the relative production of neutron stars and 
black holes. There is uncertainty in the threshold mass for 
core collapse into a black hole. Lower mass heHum star sys- 
tems should also be considered dBrown et all Il999t) . A fun- 
damental aspect in this context is whether, during the mass 
loss period, the helium core remains surrounded by a thick 
hydro gen mantle, so that common envel ope evolution pm- 
ceeds (iNelemans & van den Heuvell l200ll) . It is also possible 
that a neutron star spiralling inwards in a common envelope 
might rapidly accrete a nd eventually collapse into a black hole 
(IChevaheretal.Lll993h . On a statistical basis INTEGRAL wifl 
provide a clearer picture of the relative numbers of the various 
sub species of X-ray binary systems, and consequently a more 
balanced understanding of the rather murky final stages of bi- 
nary star evolution and by implication eventually more reliable 
estimates for the relative numbers of black hole and neutron 
stars. 
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